NGC 300 OT2008-1 AS A SCALED-DOWN VERSION OF THE 
ETA CARINAE GREAT ERUPTION 
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ABSTRACT 

We propose that the intermediate luminosity optical transient 
NGC 300 OT2008-1 was powered by a mass transfer episode from an ex- 
treme Asymptotic Giant Branch star to a Main Sequence companion. We find 
a remarkable similarity of the shapes of the light curves of the several months 
long NGC 300 OT2008-1 outburst, of the three months long 2002 enigmatic 
outburst of the B star V838 Mon, and the twenty-years long Great Eruption 
of the massive binary system Eta Carinae that occurred in the 19th century. 
Their similar decline properties hint to a common energy source: a gravitational 
energy that is released by accretion onto a main sequence star. These events 
populate a specific strip in the total energy vs. outburst duration diagram. The 
strip is located between novae and supernovae. We add recent transient events 
to that diagram and find them to occupy the same strip. This suggests that 
some intermediate luminosity optical transients are powered by accretion onto 
a compact object (not necessarily a main sequence star). These transients are 
expected to produce bipolar ejecta as a result of the geometry of the accretion 
process. 

Subject headings: stars: variables: other — stars: individual (NGC 300 OT) — 
stars: winds, outflows — stars: mass loss — supernovae: general 

1. INTRODUCTION 

The Intermediate Luminosity Optical Transient (ILOT) NGC 300 OT2008-1 (hereafter 
NGC 300 OT) was discovered in Apr 24, 2008 (Monard 2008), a few weeks after its assumed 
outburst. The distance to the host galaxy NGC 300 is ~ 1.88 Mpc (Gieren et al. 2005; 
Freedman et al. 2001). With a bolometric luminosity of Lboi = 1.6 x 10^^° erg s~^ at discovery 
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(Bond et al. 2009), it is located between supernovae (SNe) and novae, and cannot be 
easily grouped within any subtype of these two classes (Kulkarni & Kasliwal 2009). The 
pre-outburst progenitor system was discovered by Prieto (2008). It was enshrouded by 
dust (Bond et al. 2009; Berger et al. 2009;), and had a luminosity of about 6 x lO^L©, 
corresponding to a M = 10 — 15Mq star, probably in the extreme Asymptotic Giant Branch 
(AGB) stage (Thompson et al. 2008). A more massive red supergiant of mass M = 12 — 
25M0, as found by Gogarten et al. (2009) based on stellar evolution considerations, would 
also be consistent with the data. 

Bcrgcr ct al. (2009) could fit the UV-visual spectral energy distribution (SED) during 
outburst with a blackbody spectrum with a temperature of = 4670 ± 140 K, and a radius 
of Rb — 14.7 ± 1.4 AU. But if inter- and circumstellar extinction is considered (Bond et al. 
2009), based on the appearance of the absorption spectrum at maximum, a temperature of 
~ 7500 K is more consistent with the SED. The source was not detected neither in the X-ray 
band {F^ < 1.2 x 10~^^ erg s~^ cm~^) nor in the radio band. Bcrgcr ct al. (2009) attributed 
the ~ 10^ km s~^ red wing of the Ca 11 H&K absorption lines either to an infalling gas 
from a previous eruption or to a wind of a companion star. In either case the star accreting 
this matter is likely to be a main sequence (MS) star. Together with the evidence for the 
supergiant nature of the progenitor, this implies that there are two different stars in the 
system. Bond et al. (2009) interpreted the Hydrogen Balmer lines' and the Ca II IR triplet's 
double features as indicating the presence of a bipolar outflow expanding at a velocity of 
~ 75 km s~^. None of the observed hues exhibited P-Cygni profiles or velocities exceeding 
10^ km s~^ (Berger et al 2009); the latter property is a clear departure from typical SN 
behavior. 

According to Prieto et al. (2009) the ILOT event that most resembles NGC 300 OT is 
SN 2008S (Prieto et al. 2008; Wesson et al. 2009). Both ILOTs were the result of an energetic 
eruption in a dust-enshrouded 10 — 2OM0 star that survived the eruption (namely, SN 2008S 
was not a SN at all). Smith et al. (2009) suggested that the physical mechanism which 
produced SN 2008S is a super-Eddington wind, similar to the super-outbursts of massive 
Luminous Blue Variables (LBVs). Botticella et al. (2009) suggested that the progenitor 
was an extreme ( "super" ) AGB star. Bond et al. (2009) suggested that both SN 2008S and 
NGC 300 OT originated form evolved massive stars on a blue loop to warmer temperatures, 
and were subjected to increased instability due to prior mass loss. 

An asymmetric dusty environment extending a few-thousand AU surrounding NGC 300 OT 
(Patat et al. 2009) hints to a previous possible eruption. This asymmetry may further hint 
to the presence of a companion star, although we note that up to date there has been no 
definite observation proving a companion existence. Thompson et al. (2008) suggested that 
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ILOTs occur due to single star processes, e.g., electron-capture SN, an explosive birth of a 
massive WD, or an enormous outburst of a massive star. In their model for NGC 300 OT 
and SN 2008S the progenitors were luminous (~ 4 — 6 x 10^ Lq) dust-enshrouded stars, at 
the end of their AGB stage. 

Based on the model proposed by Soker (2004), we examine in this Letter (Section 
12]) whether an eruptive mass transfer onto a companion can account for the properties of 
NGC 300 OT. In that model most of the energy of the outburst was gravitational energy 
released by ~ few x O.IM© accreted by a MS B-type companion. The binary system sur- 
vives the event. In Section [3] we suggest a physical mechanism to account for most of the 
gravitationally-powered outbursts. 

2. ACCRETION BY A SURVIVING COMPANION 

We examine a model where the source of the mass is an extreme-AGB star, while the 
main source of the energy is a gravitational energy released by the mass accreted onto a MS 
companion. 

Based on previous papers (Berger et al. 2009; Bond et al. 2009; Gogarten et al. 2009; 
Prieto et al. 2009) we scale the mass of the extreme-AGB star by Mi ~ 15Mq. We also 
assume that the main sequence companion is not much lighter than the primary, and scale 
it with M2 — 8Mq, corresponding to a MS radius of R2 — 3.5Rq (our model can work for 
3Mq < M2 < IOM0). We attribute the slow outflow of ~ 50 - 75 km s"^ that probably 
contains most of the mass to the extreme-AGB star, as it is about equal to its escape speed 
(assuming a radius of i?i ~ 5 AU). The faster outflow of up to ~ 600 km s~^ fits better the 
escape velocity from the MS companion. 

Our scenario starts with some kind of instability in the extreme-AGB star that causes a 
mass of several XO.IM0 to be lost from the star at a slow velocity. In our model the instability 
is not the source of the extra energy. Therefore, we can assume that the instability does 
not increase much, or even reduces, the primary luminosity. We do not specify the source 
of the instability, but it might be, for example, a strong magnetic eruption. AGB stars are 
known to have extensive convective region with strong convection (the convective cells have 
a relative high velocity and long mixing length). If the star has a non- negligible rotation 
due to its tidal interaction with the companion, then a strong magnetic activity might be 
expected (e.g, Garcia-Segura et al. 2001). The magnetic eruption causes the primary to 
overfill its Roche lobe, resulting in both mass loss from the system and mass transfer to the 
companion. 
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According to Bond et al. (2009), the transient at maximum light exhibited a spectrum of 
an F superigant. The underlying radiation source during th.e event could be even hotter, if the 
observed spectrum was produced in an optically thick wind. However, this cannot be used to 
constrain the spectral type of the progenitor before the event. Wc take the most conservative 
approach, and take the star to have the lowest temperature possible for its mass at this 
evolutionary stage (extreme- AGB), ~ 3500 K. For an extreme AGE effective temperature 
of ~ 3500 K the radius of the progenitor of NGC 300 OT was i?i ~ 3 AU. The dynamical 
time scale for a mass of Mi — 15Mq is ~ 5 months. The outburst duration of 80 days can 
be understood as dynamical time scale. Namely, in our model a dynamical instability lead 
to high mass transfer episode. The super-Eddington luminosity of the accreting secondary 
might have helped in terminating the high mass transfer rate. For an efficient accretion the 
companion in our model of NGC 300 OT has to be very close to the primary, ~ 2Ri. We 
note that if the orbit is highly eccentric, say e = 0.9 as in rj Car, and periastron passage is 
at ttp = 2Ri 6 AU, then the orbital period is ~ 100 yr. Furthermore, if the outburst was 
caused by the periastron passage, then some high mass loss episode could have occurred 100 
years ago, but not necessarily as strong. Assuming a velocity of 75 km s~^ (as observed for 
the present bipolar outflow) these ejecta are at a distance of ~ 1500AU. If the progenitor 
was indeed an F star, its higher temperature would imply a smaller radius, and hence a 
shorter dynamical timescale of about half a month for an effective temperature of 7500 K. 
This will considerably ease the constraints on our model, as the outburst in our model is 
limited from below by the dynamical time scale. 

In our scenario most of the outburst energy comes from accretion onto the companion. 
Therefore, we can use the observed luminosity and total energy to constrain the accretion 
mass and time scale. The total radiated energy is E'rad = 10^''' erg (Bond et al. 2009; Bcrgcr 
et al. 2009). Adding the kinetic energy of the ejected mass doubles (or even more) the 
total eruption energy, which we scale with E'tot = 2 x 10*^'' erg. The accreted mass onto the 
companion is constrained to be 



where x is the efficiency of converting gravitational energy to eruptive energy. If the accreting 
star does not rotate fast, then x is very close to 1. Here we take a conservative approach, 
and scale with x = 0.5. 

The time scale of the high state of the eruption was ~ 60 days (Bond et al. 2008). This 
would give an average mass accretion rate of > O.3M0 yr~^, depending on x a-nd the other 
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parameters. 

There are some similarities between our suggested model of NGC 300 OT and the 
merger model of the outburst of V838 Mon (Tylenda & Soker 2006). In both ILOTs the 
energy is obtained from a process of mass transfer. In V838 Mon, according to the model 
proposed by Tylenda & Soker (2006), it was a merger process, in which a ~ O.3M0 star 
merged with a more massive one, on a timescale of ~ 80 days. This results in an accretion 
rate of ~ IAMq jt'^, close to that of NGC 300 OT. The ejection velocities in V838 Mon 
(~ 300 km s~^) are also similar. The duration of the V838 Mon eruption is also quite similar 
to that of NGC 300 OT. Indeed, the total energy of the two ILOTs is few xlO^^ erg. 

In some sense the NGC 300 OT eruption is a scaled-down version of the Great Eruption 
(GE) of r] Car. The GE of t] Car in 1843 ejected more than lOM© (Smith et al. 2003) from 
the 100 — I5OM0 progenitor LBV star. Soker (2004) suggested that a large fraction of this 
mass was accreted by the secondary star, of mass ~ 25Mq. Part of the mass transferred 
to the companion was expelled in a bipolar outflow, that shaped the bipolar structure of r] 
Car, the Homunculus (e.g., recent review by Smith 2009). In rj Car more data are available 
from observations, e.g. we know the binary period and the ejecta energy, therefore the 
accretion-powered model could be worked out quantitatively in more detail (Soker 2004). 

It is possible that the trigger for the impulsive mass loss episode in r] Car and NGC 300 OT 
was the same. Harpaz & Soker (2009) suggested that the mass loss from the primary in the 
GE of f] Car was triggered by magnetic eruption in the primary envelope. The magnetic 
eruption does not increase much the luminosity, but mainly causes an impulsive mass loss 
episode from the primary. The secondary accretes a large fraction of the outflowing mass. 

We now turn to discuss the light curves (Fig. [1]). The early light curve of V838 Mon 
showed 3 peaks (Munari et al. 2002), and then had a rapid decay by 5 magnitudes in about 
20 days, shorter than the ~ 80 days duration of the high-brightness phase. It is unknown 
if such peaks existed in NGC 300 OT (Bond et al. 2009), but this may be due to a poor 
time coverage of the early light curve. However, in the decline phase, the light curve of 
NGC 300 OT resembles those of V838 Mon and t] Car GE (Fig. [I]). The similarity of 
NGC 300 OT to V838 Mon and few other ILOTs was already noticed by other authors (e.g. 
Berger et al. 2009, Smith et al. 2009). 

In the case of NGC 300 OT, 120 days were required for a decay by 5 magnitudes, a time 
longer than the high-brightness phase. For the GE of r] Car, it took no less than ~ 18 yr for 
the light curve to decay by 5 magnitudes (Humphreys et al. 1999). 

Remarkably, when we normalize the timescales of the three objects without changing 
the magnitude scale we find that the light curves look very similar (Fig. [2]). A similar 
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procedure is known as "stretch correction" when apphed to SN hght curves. The normahzed 
slope of the dechne behaves the same way for all three eruptions. We consider this a hint 
that all three eruptions are governed by the same physical mechanism. For example, the 
residual energy source is not as high (even when the bolometric light curve is considered) 
as in radioactive decay in SNe. Weaker contribution is possible, e.g. from mass accretion 
at a low rate (like in the scenario for V838 Mon; Tylenda 2005). Our leading candidate for 
the powering mechanism is dissipation of accretion energy. We consider the similar slopes 
as supporting argument to NGC 300 OT being a scaled-down version of the GE of i] Car. 
The rapid decay in luminosity is determined by the decrease in the supply of accreted mass. 
Later, the slower decay is dictated by the settling of the inner part of the inflated envelope on 
the accreting object. In this later phase the luminosity is below the Eddington luminosity. 
The time scale itself depends on the accreting object and on the structure of the inflated 
envelope. 

On the other hand, when we compare NGC 300 OT to SNe (Fig. [T]) we find that the 
slope of the decline phase of NGC 300 OT is steeper than the slope of the decline phase 
of SNe. Even SN 2008S which is considered the "twin" of NGC 300 OT (Prieto et al. 
2008) shows this discrepancy - its decay was more gradual than that of NGC 300 OT in 
the V-band. This is despite the fact that its progenitor had properties similar to those of 
NGC 300 OT, and the outburst duration was similar. If SN 2008S is indeed a SN, as argued 
by Botticella et al. (2009), its kinetic energy is expected to be ~ 10 — 100 larger than the 
radiated energy. In such a case SN 2008S would be located in the region of exploding stars 
in Figure [31 consistent with its claimed SN nature. 

To summarize this section, we suggest that the eruption of NGC 300 OT can be under- 
stood as a scaled-down version of the GE of r] Car, in terms of a similar dominating physical 
process — accretion. 

3. THE TOTAL ENERGY BUDGET OF ACCRETION-POWERED ILOTS 

The objects presented in Section |2] as examples of the accretion-powered transient sce- 
nario fall between the regions of novae and SNe in the peak visible brightness vs. event 
timescale observational diagrams (e.g., Kulkarni & Kasliwal 2009; Rau et al. 2009). This 
region contains more objects that are termed ILOTs. With the accretion scenario at hand, 
it may be instructive to plot a diagram based on the total energy of the transient events 
instead of their luminosity. 

Part of the energy liberated by accretion is not radiated, but rather channeled to other 
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processes, such as lifting envelope material. Energy radiated outside the optical bands is 
usually not observed. Also, the kinetic energy of the ejecta is hard to estimate due to model- 
dependent mass loss values, and a somewhat wide range of measured velocities. Therefore, 
there is no unique way to define the total energy of the event. Nevertheless, we define the 
total 'observed' energy as the bolometric energy (estimated from the optical measurements) 
emitted during the outburst, plus the ejecta kinetic energy. The duration of the outburst is 
formally defined as the time by which the event brightness (in the V band) drops by 2mag. 
In Figure [3] this total energy is plotted as a function of the event timescale for several ILOTs 
and compared to the values for SNe and classical novae. Table [1] presents the data used to 
plot Figure [3] with references. 

SNe Ib/c, and II are grouped as Exploding Massive Stars (this includes hypernovae). 
Notice that for classical novae we plot both the observed range from Delia Valle & Livio 
(1995) (the two wavy green lines), and the model values from Yaron et al.(2005) (red dots). 
The nova model values include the kinetic energy of the ejecta, but this is usually much 
smaller or at most equal to the radiated energy (Epelstain et al. 2007). A few models from 
Yaron et al. (2005) with extreme, unobserved, values of average velocity (> 2500 km s~^) 
and mass loss (> 3 x 1O~^M0) were not included in the plot. 

It is evident that the locus of ILOTs in Figure [3] is more separated from that of classical 
novae than in the luminosity vs. timescale plane (e.g., Kulkarni & Kasliwal 2009; Rau et al. 
2009). The difference from SNe is also clearer, even though less striking. The use of energy 
instead of luminosity apparently highlights the physical mechanisms underlying the various 
transient classes. The region of the diagram populated by the supposedly accretion-powered 
events is tentatively marked with a colored strip. We take the appearance of this sequence 
as further evidence for a common physical mechanism behind these transients (although not 
all ILOTs are necessarily powered by accretion). As for P Cygni, which is not yet proven 
to be a binary system, we predict that a binary companion exists in that system, or that a 
common envelope occurred in the recent past.). 

To summarize, we suggest that most (but probably not all) of the ILOTs are accretion- 
powered events. Based on the common energy source they should be grouped in one class 
together with stellar merger tidal, disruption flares, and large eruptions in LBV binary 
systems. The shape of the light curve (Fig. [T]), related to the accretion as the energy source, 
is an important distinguishing feature of many systems in this group. 

We also predict that because of the geometry of the accretion process the ejecta in this 
class of transients should systematically exhibit bipolar structure, as suggested by Soker 
(2004). We note that in some single star models circumstellar material can also posses 
bipolar structure (e.g.. Smith 2009). 
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Table 1: Data used to plot the total observed energy vs. event timescale in Figure [31 Total 
energy included radiative and kinetic energy. References: 1. Berger et al. (2009), 2. Tylenda 
& Soker (2006), 3. Kulkarni et al. (2007a), 4. Kulkarni et al. (2007b), 5. Rau et al. (2007), 
6. Pastorello et al. (2007), 7. Ofek et al (2007), 8. Smith et al. (2009), 9. Soker, Frankowski 
& Kashi (2009), 10. Barbary et al. (2009) 11. Rich et al. (1989), 12. Mould et al. (1990), 
13. Foley et al. (2009), 14. Valenti et al. (2009), 15. Arnett & Fu (1989), 16. Bethe & 
Pizzochero (1990), 17. Davidson & Humphreys (1997), 18. Smith et al. (2003), 19. Frew 
(2004), 20. Humphreys et al. (1999), 21. Smith (2005), 22. Richardson et al. (2009), 23. 
Smith & Hartigan (2006), 24. Rau et al. (2009), 25. Branch (1992), 26. Richardson et al. 
(2006), 27. Yaron et al. (2005). 



Name 


Dm'ation 
[d] 


Total energy 
[10^8 erg] 


rcf. 


Remarks 


NGC 300 OT 


80 


0.2-0.4 


(1) 




V838 Mon 


90 


0.3-1 


(2) 




M85 2006 OT (SN) 


180 


50-100 


(3,4,5,6) 


Assuming it is a weak SN 


M85 2006 OT (Acc.) 


180 


0.1-0.8 


(6,7) 


Assuming it is accretion-powered 


SN 2008S 


75 


~ 1.3 


(8) 


We added kinetic energy equal to the 
radiated energy (upper estimate from 
(7)) 


SCP 06F6 


150 


- 50 


(9,10) 


Assuming the source is extragalactic; 
we include the rising phase because the 
profile is symmetric 


M31 RV 


70 


0.9 


(11,12) 




SN 2008HA 


32 


2-50 


(13,14) 




SN 1987A 


130 


1500 


(15,16) 




r] Car GE 


7300 


80 


(17,18) 




7? Car LE 


1800 


3 


(19,20,21,22) 


Lesser Eruption of i] Car (1890) 


P Cygni 1600 AD 


2200 


2.5 


(19,20,23) 




SNe la 


30-70 


500-2000 


(24,25) 


Radiated energy is negligible 


Exploding Massive Stars 


20-400 


100-10^ 


(24,26) 


SNe Ib/c, and 11 including hypernovae 


Classical Novae 


1-1000 


10-^^-10-2 


(3,27) 


Dots are models from (24), the wavy 



band is computed using luminosity and 
duration from (3) 
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Fig. 1. — Comparison of the V-band light curves of the 77 Car Great Eruption (GE), V838 Mon and 
NGC 300 OT. The timescale was norniahzed so that 1 time unit equals 1 year for 77 Car GE, 2.2 days for 
V838 Mon, and 5.6 days for NGC 300 OT. For NGC 300 OT the R-band is also plotted, for which there 
has been one observation before the maximum (Bond et al. 2009), marked with a red circle. Top: the three 
separated light curves; the apparent V-mag axis was not rescaled. Bottom: The same curves translated 
vertically to bring peak luminosities to overlap (see legend for the shift values). It can be easily seen that 
the slope of the decline phase and its rate of change are similar for the three eruptions. 



SNII 

SN2008HA 

SN la 

SN 2008S 

NGC 300 OT 

-10 10 20 30 40 50 60 70 

time 

Fig. 2. — Comparison of the light curve of NGC 300 OT to SN light curves. The timescale was normalized 
so that 1 time unit equals 5.6 days for NGC 300 OT, 4.5 days for SN 2008S, 1.3 days for SN 2008HA, 2.9 
days for SN la, and 1.3 days for SN II. The magnitude was shifted to bring peak luminosities to overlap. It 
can be easily seen that the decline of NGC 300 OT is much steeper than the decline of SNe. The references 
for the data are: SN 2008S - Botticella et al. (2009); SN la (represented by SN 1991T) - Patat et al. (2001); 
SN II (represented by SN 1993J) - Patat et al. (2001). 
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Fig. 3. — The total (radiated plus kinetic) energy of the transients discussed in Section [3] as a function of 
the duration of their outbursts. We did not include the energy which is deposited in lifting the envelope since 
it is not observed. See text and Table 1 for information about the different transients and the derivation of 
their total energy. The shaded area marks the group of transients we study in the paper. We propose that 
most of these transients are associated with mass-transfer processes. 



